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Abstract. In this paper we give a conjectural refinement of the Davcnport- 
Heilbronn theorem on the density of cubic field discriminants. We explain how 
this refinement is plausible theoretically and agrees very well with computa- 
tional data. 



1. Introduction 

Let a„ be the number of isomorphism classes of abehan cubic fields with discrim- 
inant n. Let bn be the number of isomorphism classes of non-abelian cubic fields 
with discriminant n. The numbers a„ are very well understood. The numbers 6„ 
have been the subject of extensive theoretical and computational study for at least 
sixty years, but are less well understood. The object of this note is to contribute 
to the study of these 6„, by bringing together the theoretical and computational 
literature. 

For a E {— , +} , define 

X 

(1.1) ga{x) ^y^^bgn. 

n=l 

Put C- =3/4 and C+ = 1/4. The main theoretical result is due to Davenport and 
Heilbronn [8]: 

(1-2) ga{x) - C'"^^^;. 

The numbers ban have been computed for larger and larger n by many authors. 
The most recent results are due to Belabas [3], who introduced a vastly improved 
method and worked through the cutoff x = 10^^. 

It has been discussed by several authors that the "match" between computational 
evaluation of ga{x) and the asymptotic formula (1.2) is not very good, as illustrated 
by Tables 1, 2, and 3 below. For example, the paper [9] dealt with the case a = —; 
here is an excerpt from pages 322-324, with some trivial notational changes: 

If, however, the reader were to plot the density versus x, he would be 
somewhat astonished to see that this density is increasing so slowly 
that his first impression would be that it will not make it to Davenport- 
Heilbronn limit. Thus it remains a challenging problem, assuming that 
the D-H limit is not in error, to explain the origin of this very slow 
convergence. This problem was indicated by Shanks in [13], [1], and 
on the real side in [13], [10], where the problem is further aggravated 
by even slower convergence. To date, and to our knowledge, no good 
quantitative explanation of this phenomenon has been given. 

In this paper we conjecture a quantitative explanation of this slow-convergence 
phenomenon. 
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2. The abelian case 
It is best to briefly review the situation in the abelian case. One has the formula 



(2.1) 



Put 

X 

(2.2) /(a;) = ^a„. 

n=l 

Cohn [4] showed that the function (2.1) continues analytically to Rc(s) > 1/4, 
except for a simple pole at s = 1/2. He computed the residue of this pole. This 
gave him more than enough information to conclude that 



,2.3) m-l^f n 



p=l (6) 



Define h-{x) = g-{x) and 
(2.4) h+{x)=g+{x) + ^f{x). 

From (1.2) and (2.3) one sees that (1.2) would remain true with replaced by /iq.. 
We henceforth focus more on ha, as it appears naturally in (3.2) and (3.3). 

3. The basic conjecture 

Belabas proved [3, (8) and (20)] that, for either a, the difference between the 
sides of (1.2) is o(a;/(loga;)^). Prom his extensive computational data he speculated 
[3, page 1232] that these differences arc o(.x/(log.T)'') for arbitrarily large (3. We 
believe a stronger statement. Put K- = 3/(3 + v^) and Kj^ = •\/3/(3 + v^). 
Conjecture. For each a, 

(3 1) h(x)-C ^x + K C(l/3)3(3 + 73)r(l/3)^ 3/6 + 

(3.1) ha[x)-Ca^^^^^x + K„ 107r3C(2)C(5/3) ^ "^'^ 

The heuristic we arc about to sketch is too crude to suggest any particular error 
term. The o{x^/^) we have written corresponds to the weakest conjecture in which 
our conjectured new term actually plays a role. However, Tables 1 and 2 below 
suggest that the error term should be reducible at least most of the way to o{x^/'^). 
We are led to our conjecture theoretically as follows. Define a Dirichlet series 

Here is the set of isomorphism classes of cubic fields k with discriminant of sign 

a. Also, Dk is the discriminant of k while Aut(fc) is the automorphism group of 
fc; so I Aut(fc)| is 3 or 1 according to whether k is abelian or not, the former case 
occuring only for a = +. 

For S a finite set of primes, define a weighted version of (3.2) 

Here the weights rjk^s{s) are given on page 122 of [6]. These weights are themselves 
Dirichlet series ^ Cm{k, S)m~^ with integral coefficients Cm{k, S) decreasing to the 
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Kronecker symbol 6mi as S grows. Thus ^a,s in the sense of formal Dirichlet 
series. 

The point of introducing the approximations £,a,S is that one understands them 
analytically better than one understands itself. From Davenport's earlier work 
[7] one knows that the rightmost pole of ^q, j is at s = 1. The analytic behavior of 
^a,0 was investigated more explicitly by Shintani [14], who found the next pole at 
s = 5/6. More generally, it is proved in [5, Prop. 6.2], that each ^„,s continues to 
the whole ,s-planc. mcromorphically with only poles at s = 1 and s = 5/6. These 
poles are simple, their residues Caries and K^r^/Q^s being deducible from the first 
two equations in [5, Thm. 6.2i]. In the case s = 1, one has to set 6^ equal to zero 
in the corresponding equation; this h„ corresponds to quadratic fields, included in 
£,a,s in [5] and [6], but not included here. 

Prom the formulas in [5, Thm. 6.2i], one sees that the numbers ri,5 are positive, 
and decreasing as S grows. Similarly, the numbers r^^/Q s arc negative and increas- 
ing as S grows. Let ri and r^/Q be their respective limits. One does not have enough 
control at present to rigorously apply any "poles-control-growth" theorems. How- 
ever the general formalism suggests that ha{x) grows as CaVix + Ka{'o/^)rr-^iQX^f^ . 
This is the source of our Conjecture 3.1. The constant rs/g is actually nicer than we 
had initially expected from its source as the limit of the r5/6,s; special zeta values 
arise because of the factorization of the right side of (5.3) below. 

To follow our heuristic in a more detailed way, one should first understand it 
in the more refined context of Section 5 below. We recommend looking next at 
Prop. 2.1 and the discussion aroimd it in [6]; this treats in more classical language 
the special case of [5, Thm 6.2i] we need, except it does not repeat the formula for 

^5/6,S- 

4. Numerical evidence for the basic conjecture 

Our conjecture for a = — compares with Belabas's data [3, page 1232] for 
negative-discriminant cubic fields as in Table 1. Here and below, x = 10^. Also 

Table 1. Computation and theory: negative-discriminant cubic fields 



j 


h-{x) 


h_ix) 


/i_(x) 


h-ix)-Hl{x) 


H^{x) 


H*_{x) 


2.1/2 


2 


7 


0.337 


0.7843510 


-0.1925 


3 


127 


0.611 


0.9993210 


-0.0027 


4 


1,520 


0.731 


0.9943300 


-0.0867 


5 


17,041 


0.819 


0.9998781 


-0.0066 


6 


182,417 


0.877 


1.0001096 


0.0200 


7 


1,905,514 


0.916 


1.0000100 


0.0060 


8 


19,609,185 


0.943 


0.9999394 


-0.1188 


9 


199,884,780 


0.961 


0.9999850 


-0.0951 


10 


2,024,660,098 


0.974 


1.0000009 


0.0176 


11 


20,422,230,540 


0.982 


1.0000003 


0.0218 



Ha{x) denotes the right side of (1.2) while H^{x) denotes the first two terms of the 
right side of (3.1). Clearly Hl_ matches the data h- substantially better than iJ_ 
does. The poor match of H- to h- in the range 10^ < a; < 10^ formed the basis of 
the quote in Section 1. 

For a = -h, we start from the values of g+{x) + f{x) listed on [3, page 1232]. We 
use (2.1) to compute f{x) and thereby deduce g+{x). Our hope that the error term 
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Table 2. Computation and theory: positive-discriminant cubic fields 



J 


1 (tA 


fix) 


H+{x) 


h+{x) 
H*+{x) 


h+ix)-Hlix) 

^1/2 


g+ix)-Hl{x) 
a;i/2 


2 





2 


0.096 


8.5889786 


0.0589 


-0.0078 


3 


22 


5 


0.341 


1.0461129 


0.0330 


-0.0197 


4 


366 


16 


0.536 


0.9900166 


-0.0374 


-0.0908 


5 


4,753 


51 


0.688 


1.0010833 


0.0163 


-0.0374 


6 


54,441 


159 


0.786 


0.9988436 


-0.0631 


-0.1161 


7 


592,421 


501 


0.855 


0.9999134 


-0.0162 


-0.0690 


8 


6,246,698 


1,592 


0.901 


1.0000259 


0.0161 


-0.0369 


9 


64,654,353 


5,008 


0.933 


1.0000097 


0.0198 


-0.0330 


10 


661,432,230 


15,851 


0.954 


0.9999988 


-0.0081 


-0.0609 


11 


6,715,773,873 


50,152 


0.969 


1.0000002 


0.0046 


-0.0482 



in Conjecture 3.1 can be reduced to o{x^^'^) would require that the numbers in the 
last column of Table 1 and the next to last column of Table 2 tend to zero. The 

data seems rather borderline in support of o{x^^'^), but strongly indicative that our 
conjecture as stated is true, and moreover that the o{x^^^) can be reduced at least 
somewhat. 

From Cohn's result (2.3), at most one linear combination g+{x) + uf{x) can 
agree with H^{x) to within o(a;^/^). The last two columns of Table 2 suggest that 
u = 1/3 is more likely than u = 0. In fact u = 0.31 . . . would give the best fit, in a 
least squares sense, and this agrees well with the theoretically favored u = 1/3. 



5. The refined conjecture 

The Davenport-Heilbronn theorem comes in a more refined version and so does 

our conjecture. Let A be the set of symbols {111, 21, 3, 1^1, 1^}. Then for K a 
cubic field and p a prime, the splitting behavior of p in K determines an element 
of A. Up is unramified, then it either splits completely (111), splits partially (21), 
or remains inert (3). If p is ramified, then it is either partially ramified (1^1) or 
totally ramified (1^). 

Up until now, we have been using a to denote a sign, -|- or — . Henceforth, we use 
a to denote the following data. First, a finite set support(Q!) C {oo, 2, 3, 5, ... } of 
places of Q. Second, if oo e support(a), a sign aoo G {+, — }, as before. Third, for 
each prime p G support (a), an element ap G A. Let be the set of isomorphism 
classes of cubic fields meeting the given local condition a„ for each v G support (a). 
One naturally has fa, gai ha, Ca, and £,a,S, all direct generalizations of the notions 
introduced above in the special case support (a) = {oo}. 

Define Coc.- = 3/4 and Coo,+ = 1/4, as in Section 1, and Koo,- = 3/(3 -|- V^) 
and K^^+ = V3/(3-F V3), as in Section 3. Define p-adic analogs as follows: 



(5.1) 



ap 




K 


111 

21 
3 

1^1 

13 


l/QCp 
l/2Cp 

i/iCp 

l/pCp 
l/p^Cp 


(i+p-y^f/QKp 

ii+p-'/^)ii+p-y^)/2Kp 

{l+p-')/3Kp 

{i+p-'/^r/pKp 

(l+p-V3)/p2^^ 
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Here 
(5.2) 

(5.3) 



(l_p-5/3)(i+p-l) 



l-p- 



-1/3 



axe normalizing factors, chosen so that the column-sums ^ Cp,ap and J2 ^p,ap are 
both 1. The values of Cp^ap are taken from e.g. [10, page 589] while the Kp^ap are 
taken from [5, Table 5.1]. Note that both the Cp_ap and the Kp^^p are actually 
given by procedures uniform in ap. For example, the ratios [Kp^apKp) / {Cp^apCp) 
are all given by replacing a factor in Up by a factor (1 +p^^^), e.g. 1^1 ^ 

(l+pV3)(l+pl/3). 

We expect that Conjecture 3.1 holds verbatim for these more general a, with 

tJ€:Support(Q) 



(5.4) 

(5.5) = \[ K^,ap. 

i'Gsupport(a) 

Our heuristic continues to make sense in this more localized context. In fact the 
basic references [5] and [6] are in this context. 

6. Numerical evidence for the refined conjecture 

In this section, wc work with twenty different a. All of them have support {oo,p} 
and = +• We let p run over {2, 3, 5, 7} and let ap run over the five-element set 
A. For simplicity, we compare computation with theory only at a; = lO''. 

Wc get the numbers ga(lO^) by rcnormalizing the first four lines of [10, Table 11]. 
We get the numbers /^(lO'') by a simple procedure based on (2.1). Including these 
extra 501 fields, with total mass 501/3 = 167, improves the fit slightly, as would be 
expected from line j = 7 of Table 2. 

Table 3 gives the comparison. This table illustrates how the slowness of the 

Table 3. Computation and theory: j>adic decomposition (a; = lO'') 





p= 


2 


P= 


3 


P= 


5 


P= 


7 




ha{x) 


ha{x) 


ha{x) 


ha{x) 


ha{x) 


ha{x) 


ha (x) 


ha{x) 






H*{x) 


Haix) 


H*{x) 


Haix) 


H*(x) 


Hjx) 


H*{x) 


111 


0.706 


0.9999 


0.721 


0.9998 


0.734 


0.9994 


0.740 


0.9994 


21 


0.851 


0.9999 


0.856 


1.0001 


0.858 


0.9996 


0.858 


0.9997 


3 


0.924 


1.0001 


0.923 


0.9999 


0.920 


1.0002 


0.917 


1.0008 


IH 


0.836 


0.9996 


0.835 


0.9997 


0.832 


1.0002 


0.829 


0.9985 


13 


0.909 


1.0003 


0.903 


0.9998 


0.895 


1.0006 


0.891 


1.0025 



convergence ha{x) / Ha{x) to 1 seems to depend systematically on ap, the totally 
split case ap = 111 being the slowest. On the other hand, ha{x) / H*{x) is much 
closer to 1, with no apparent tendency to be above or below. 

7. Three concluding remarks 

By using the five clcmcmt set A, wc arc falling slightly short of the most locally 
refined conjecture, because there are different ways a prime can be ramified. Thus, 
for example, if a field K isl^ at p = 3, then there are actually nine possibilities for 
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its 3-adic completion K3; both C^^is and -fl'3,13 should be split into nine parts, in 
the proportions indicated in the {p, Up) = (3, 1^) slot on the following chart: 





ap 


p = 2 


p = 3 p > 5 


(7.1) 


1^1 


111111 

4 4 8 8 8 8 


11 11 
22 22 




13 


1 


111111111 1111 

3 3 9 81 81 81 81 81 81 3 3 3 



Similarly, if a field K is 1'^ at p > 5, there are either 1 or 3 possibilities for K3, 
according to whether p = 5, or 1 modulo 6. In general, a p-adie algebra Kp of type 
ap contributes to Cp,ap and Kp^ap in proportion to 1/(| Aut{Kp)\\DKp\)- 

Our basic theoretical references [5] and [6] work over arbitrary number fields as 
well as over finite extensions of Fp{t) for p > 5. We expect that the analog of our 
fully refined Conjecture 3.1 holds in this extra generality too, but at present there 
is no corresponding published computational data. 

The most important remaining problem is to prove our conjecture. The generally 
pessimistic discussion on [6, page 124] and [2, page 620] suggests to us that the 
way may be difficult. However, one ingredient of a proof might be the functional 
equation of ^a,s with respect to s 1-^ 1 — s, studied in [14], [5], and [11]. Another 
ingredient might be [12, Thm. 3], which concerns growth of arithmetic functions 
whose associated Dirichlet series satisfy such a functional equation. 
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